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Abstract 

This paper studies the possibility of upper bounding the position error of an estimate for range based 
positioning algorithms in wireless sensor networks. In this study, we argue that in certain situations when 
the measured distances between sensor nodes are positively biased, e.g., in non-line-of-sight conditions, 
the target node is confined to a closed bounded convex set (a feasible set) which can be derived from the 
measurements. Then, we formulate two classes of geometric upper bounds with respect to the feasible 
set. If an estimate is available, either feasible or infeasible, the worst-case position error can be defined 
as the maximum distance between the estimate and any point in the feasible set (the first bound). 
Alternatively, if an estimate given by a positioning algorithm is always feasible, we propose to get the 
maximum length of the feasible set as the worst-case position error (the second bound). These bounds 
are formulated as nonconvex optimization problems. To progress, we relax the nonconvex problems and 
obtain convex problems, which can be efficiently solved. Simulation results indicate that the proposed 
bounds are reasonably tight in many situations. 
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I. Introduction 

Recent advances in technology have instigated the use of tiny devices as sensors in large distributed 
wireless sensor networks (WSNs). A sensor device is capable to sense its environment for monitoring, 
controlling, or tracking purposes for both civil and military applications (T). Due to drawbacks in 
using GPS for WSNs, extracting the position information from the network, also called localization, 
has been extensively studied in the literature It is commonly assumed that there are a number 

of fixed reference sensors, also called anchors, whose positions are a priori known, e.g., by using GPS 
receivers Q. To find the position of other sensor nodes at unknown positions, henceforth called target 
nodes, it is assumed that there are some types of measurements, e.g., time-of-arrival, angle-of-arrival, or 
received signal strength, taken between sensor nodes (H. 

During the last decades, various positioning algorithms have been proposed in the literature. Different 
positioning approaches can be categorized based on various factors {8|. For instance, as long as an accurate 
model of measurements and the statistics of the measurement errors are known, classic estimators, e.g., 
the maximum likelihood (ML) and the least squares (LS) approaches, can be employed successfully 
to solve the positioning problem. When the distribution of the measurement errors is unknown or the 
computational complexity of classic estimators is too high, a number of simple techniques can be applied 
to the problem. For example, based on a geometric interpretation, the authors of 0, iflOl formulated 
the positioning problem as a convex feasibility problem (CFP) and applied the well-known orthogonal 
projection onto convex sets (POCS) approach to solve the problem. This method turns out to be robust 
against non-line-of-sight (NLOS) conditions ifTTTl . POCS was previously studied for the CFP and has 
found applications in several research fields lfl2ll . lfT3l . 

Positioning algorithms can be evaluated based on different performance metrics such as complexity, 
accuracy, and coverage [8). In the literature one way to assess the positioning algorithms is to evaluate 
the position error, defined as the Euclidian norm of the difference between the position estimate and the 
true position. There are a number of techniques to evaluate the performance of an algorithm based on 
the position error. For instance, a lower bound on the mean square position error is a common metric. 
There exist a number of such lower bounds in the literature, e.g., the Cramer-Rao lower bound (CRLB), 



An upper bound on position error 
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Fig. 1. An example of the application of an upper bound on the position error for traffic safety. A solid circle defines the area 
in which a vehicle definitely lies. In this figure based on an upper bound on the position error, car 2 and 3 might collide. 

which can serve as benchmarks. The CRLB, which gives a lower bound on the variance of any unbiased 
estimator, can be computed if the probability density function (PDF) of the measurement error is known 
and satisfies some regularity conditions Ifl4l . Generally, different benchmarks in the literature are used 
to statistically assess a positioning algorithm, which implies that the error in a single position estimate 
cannot be characterized in a deterministic fashion. 

Besides a lower bound on the position error, in some applications it may be useful to know the 
worst-case behavior of the position error. Such knowledge may be useful not only for evaluation of 
different services provided by WSNs but also for design and resource management 12, |fl31 . Similarly in 
evaluation of the worst-case position error, we may be interested in assessing a single point estimate. As 
an example consider Fig.Q] which shows how a nontrivial (i.e., finite) upper bound on position error can 
be used by a traffic safety application. If an estimate of a vehicle and a nontrivial upper bound on the 
position error are available, we can define an area in which the vehicle is certainly located, e.g., a disc 
centered at the position estimate and with a radius equal to the upper bound on the position error. By this 
approach, we may be able to decrease the number of collisions between vehicles. In general, computing 
the maximum possible position error might be difficult, but one may be able to derive an upper bound 
on the maximum possible position error. To the best of our knowledge, there is no specific work in the 
literature on deriving an upper bound on the position error. In this study, we aim at tackling this subject. 



In general, the concept of an upper bound on the position error (or any estimation error) seems to 
be shaky. In fact, it is not clear that is meaningful to study upper bounds, since the position error can, 
in general, be arbitrarily large. In this study, however, we argue that in some situations it is possible 
to reasonably determine the worst-case position error. For instance, if a target node position belongs 
to a closed bounded set (a feasible set), the worst-case position error can be defined with respect to 
the feasible set. For example, for distance-based positioning, if measurement errors are assumed to be 
positive, a convex set including the target node can be defined from measurements. The feasible set, in 
which the target node is located, is the intersection of a number of balls (in a 3-dimensional network) 
or discs (in a 2-dimensional network) centered at the position of reference nodes |[T6l . The assumption 
of positively biased measurement errors is fulfilled in some scenarios. For instance, in NLOS conditions, 
the measured distances are often much larger than the actual distances. Assuming a closed bounded 
(compact) convex set derived from positively biased distance measurements, a position estimate given by 
an algorithm can be either feasible or infeasible with respect to the feasible set. If an estimate is available 
(feasible or infeasible), it is reasonable to define the maximum distance from the estimate to any point in 
the feasible region as the worst-case position error. This idea yields an upper bound on the position error 
as the solution of a nonconvex optimization problem. Alternatively, a number of positioning algorithms, 
e.g., POCS, give one feasible point as an estimate. In this type of estimators, we can upper bound the 
position error as the maximum lengthj of the feasible set. To find the maximum length of the feasible 
region, we consider an outer-approximation of the feasible set and find the minimum Euclidean ball or the 
minimum ^ ball (minimum bounding box) covering the set. We further relax the nonconvex optimization 
problem and derive a convex optimization problem. Obviously, if a feasible point is available, the first 
upper bound, i.e., the maximum distance from the estimate to any point in the feasible region, gives a 
tighter upper bound compared to the second bound, i.e., the maximum length of the feasible region. 

Note that the technique introduced in this paper can be applied to every estimation problem when the 
unknown parameter vector belongs to a compact, finite- volume, convex set. 

'By the maximum length of a set, we mean the maximum £2 norm of the difference between two points (not necessarily a 
unique pair of points) in the set. 
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In summary, the main contributions of this study are: 

• introducing the concept of an instantaneous upper bound for a single point position estimate when 
the distance measurements are positively biased, e.g., in NLOS conditions; 

• proposing an upper bound on the position error based on a convex relaxation technique when an 
estimate of the target position is available (feasible or infeasible); 

• proposing three upper bounds for an estimator always giving a feasible point as an estimate (e.g., the 
POCS estimate) based on the idea of the maximum length of the feasible set or a relaxed feasible 
set including the target node. 

The remainder of the paper is organized as follows. Some preliminary requirements are studied in 
Section |II1 Section [III] explains the signal model considered in this paper. In Section ITVl a geometric 
positioning algorithm (POCS) is briefly studied. Two types of upper bounds are derived in Section [V] 
Simulation results are discussed in Section IVj Finally, Section I VII I makes come concluding remarks. 

II. PRELIMINARIES 

A. Notation 

The following notations are used in this study. Lowercase and bold lowercase letters denote scalar 
values and vectors, respectively. Matrices are written using bold uppercase letters. By nxn we denote 
the n by n zero matrix, and we use 0„ as the n-vector of n zeros. l n and I„ denote the vector of n ones 
and the n by n identity matrix, respectively. The operator tr(-) is used to denote the trace of a square 
matrix. The i v norm is denoted by || • || p . Given two matrices A and B, A >- (^)B means that A — B 
is positive (semi)defmite. S n , W 1 , and M™ denote the set of all n x n symmetric matrices, the set of all 
nxl vectors with real values, and the set of all n x 1 vectors with nonnegative real values, respectively. 

B. Quadratic ally constrained quadratic programming 

Let us consider a quadratically constrained quadratic program (QCQP) as 

maximize x T Anx + 2b^x + en 

subject to x T AjX + 2bfx + a < 0, i = l,...,N, (1) 



for Ai e b; G 1", and a € R. 

The QCQP problem £[]), in general, is nonconvex and difficult to solve except in some specific 
cases lfT7l . For the nonconvex case, there are a number of techniques to approximately solve the problem. 
One powerful approach is the semidefmite relaxation technique lTT8l - |[23l . Considering a property of the 
trace operator, i.e., x T AjX = tr(AjXx T ), the QCQP problem in £T|) can be written as 

maximize tr ( Bn \x T l] T |x T l] ) 

xGR" V L J L -V 

subject to tr (b { [x t l] T [x T l] ) < 0, i = 1, . . . , N, (2) 



where 



Ai hi 



(3) 



Now, by replacing Z = [x T l] [x T l] and noting that Z is a rank-1 symmetric positive semidefmite 
matrix, we get an equivalent problem of © as 

maximize tr (BnZ) 

subject to tr(BjZ)<0, i = l,...,N, 

Z^O, Z(n+ l,n+ 1) = 1, rank(Z) = 1. (4) 

Due to the nonconvex constraint rank(Z) = 1, the optimization problem in (0]) is still nonconvex. To 
change it to a convex problem, we drop the rank-1 constraint and obtain a semidefmite programming 
problem (SDP) as follows: 

maximize tr (BnZ) 

ze§"+! 

subject to tr (BiZ) < 0, i = l,...,N, 

Z>r0, Z(n + l,n + l) = 1. (5) 

To refer to the QPCP formulated in dH) throughout this paper, we use QP{Aj, bj, C{}^L . Similarly, to 
refer to the SDP relaxation derived in © originated from QCQP in (Q]), we use SDP{A.;, bj, C{}^L . For 
the optimal values of the objective function of the QCQP and the corresponding SDP relaxation in (Q]) 



and in ©, we use v qp {Aj, bj, Ci}f =Q and v S( jp{Aj, bj, Ci}f =Q , respectively. By adopting the relaxation, 
i.e., dropping the rank-1 constraint, we expand the feasible set, therefore, the objective function in ® is 
maximized over a larger set than in (Q]), thus 

v qp {Ai, bj, Ci}f =0 < v sdp {Ai, bi,a}f =0 . (6) 

If the rank of matrix Z for the optimal solution in (f5]) is one, then, the solution in (f5]) is equal to the 
optimal solution in £[]). In general, the optimal solution in (f5]) has rank higher than one, and then a rank- 
1 approximation can be applied to the optimal solution in Q, e.g., using a method based on singular 
value decomposition or an approach based on randomization ll20l . For details of rank-1 approximation 
techniques from a higher rank matrix, see, e.g., EUI . |[23l . ll24l . 

Note that using the Lagrange dual approach, a similar problem as the SDP relaxation in ([5]) can be 
obtained [18]. We complete this section by a simple and useful property of the quadratic inequality. 

Lemma 2.1: For a quadratic function x r Ax + 2b T x + c, where A £ b G W 1 , and c 6 K , the 
following statement always holds true: 



x T Ax + 2b T x + c > 0, Vx € 1" 

Proof: See HI. 



A b 

b T c 



y o. (7) 



C. Bounds on estimation errors given a realization of the measurement vector 

Consider an unknown parameter vector x £ W 1 . Regardless if we model x as random or unknown 
deterministic, we can define the set of the possible values of x as 

X = {possible values of x} C M. n 

Suppose m is the observed realization of the (random) measurement vector M. Given the event M = m, 
the set of possible values of x changes to 

X(m) = {possible values of x : M = m} C X. 

The estimate of x, denoted by x(m, f) € R n , is a function of the observed data m and some algorithm 
tuning parameters, e.g., initialization, step size, termination criterion, etc., which are collected in the 
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vector f. The f -vector is chosen, possibly randomly, from the set F. In other words, f G F completely 
determines how the estimator maps the observed data m to the estimate x, and the set F defines a class 
of estimators. We can now define the set of possible values of x(m, f ) when f can take on any value in 
F as 

X(m) = {x(m,f) : f G F} C R n . 
We can define three upper bounds on the £2 norm of estimation error e=||x(m,f)— x||2as 

e < tti(x(m, f)) = sup ||x(m, f) — x||2, (8) 

x6AT(m) 

e < u 2 (x) = sup ||x-x|| 2 , (9) 
x6^(m) 

e < U3 = sup ||x — x||2- (10) 

x6A"(m), xeA'(m) 

We note that all bounds depends on m, which, for simplicity, is neglected in the notation. Moreover, it is 
easy to see that iii(x(m, f)) < 113 and U2(x) < U3. Fig.|2] graphically shows the different upper bounds. 

Remark 1: the bound iii(x(m, f)) is an upper bound of the norm of the estimation error for a certain 
estimate (f and m are fixed). Hence, if ui(x(m, f)) can be computed together with the estimate, this 
would greatly increase the value of the estimate, since we can now guarantee that the norm of the 
estimation error in x(m, f) does not exceed ui(x(m, f)). This is a much stronger statement than to 
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provide a statistical quality measure, such as the mean-squared error of the estimator, 

E M {||x(m,f)-x||l}, 

where Em denotes expectation over the distribution of M. 

Remark 2: the bound U3 could potentially be computed together with the estimate and is therefore 
of value in a practical situation. However, 113 will only be interesting if it is easier to compute than 
ui(x(m, f)), since ui(x(m, f)) < U3. 

Remark 3: the bound U2(x) can be interpreted as the error of the worst estimate that is computed 
from the observed data m by the class of estimators defined by T. This is useful to judge the worst 
case performance of a class of estimators. However, since the bound is a function of x (the unknown 
parameter), it cannot be computed together with an estimate, and its practical value is therefore limited. 
We can also formulate lower bounds by replacing sup with inf in Eqs. ([8])- (fTOb . 

e >€i(x(m,f)) = inf ||x(m, f) - x|| 2 , (11) 

x€X(m) 

e > ^2( x ) — inf || x — x ||2> (12) 

xei'(m) 

e>£ 3 = inf ||x-x|| 2 . (13) 

xeA"(m), xGA'(m) 

In general, there are no guarantees that any of the bounds in Eqs. (f8T> — (fT3T> are nontrivial, i.e., that the 
upper bounds are finite and the lower bounds are greater than zero. For example, if the set X(m) or 
Af(m) is unbounded, it is clear that the upper bound ([8]) or (fTOt is trivial. However, as we will see in 
the remainder of this paper, there are indeed practical situations when the bounds are nontrivial. 

III. System Model 

Let us consider an n-dimensional network, n = 2 or 3, with N reference nodes at known positions 
a« = [tii,i • • • cLi,n] T £ I^ n > i = 1, —,N. Suppose that a target node is placed at an unknown position 
x = [xi ■ ■ ■ x n ] T E W 1 . The range measurement between the target and reference node i is given by 

di = di(x,aj) + a, i = l,...,N, (14) 
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where dj(x, a*) is the actual Euclidian distance between the target node and reference node i, i.e., 
dj(x, aj) = ||a, — x||2, and is the measurement error. 

In the literature the measurement error is commonly modeled as a zero mean Gaussian random 
variable [H, |4), ll25l . In some scenarios, however, other distributions seem to be more reasonable. 
For instance, in NLOS conditions the measured distances are larger than the actual distances with 
high probability. A number of distributions have been considered to model NLOS conditions, e.g., an 
exponential distribution or a uniform distribution ll26l . The Gaussian distribution with large positive 
mean has also been considered to model the NLOS condition |[26l . lETl . In this paper for the purpose 
of deriving the upper bound, we assume that the distance measurements are positively biased, meaning 
the measurement errors are nonnegative. The positive measurement assumption can be fulfilled, e.g., in 
NLOS conditions (with high probability). 

The positioning problem, then, is to find the position of the target node based on the positions of iV 
reference nodes and measurements made in (fl4l ). 

IV. Positioning algorithms 

A classic method to solve the problem of positioning based on measurements taken in (fl4l) is to employ 
an ML estimator if the distribution of the measurement error ej is known. Otherwise, when the statistics 
of measurement errors are unknown, one can apply the LS minimization as lfl4ll . |[28l 



The solution to (031 ) coincides with the ML estimate if the measurement errors are zero mean, 
independent and identically distributed Gaussian random variables fl4ll . In general, the LS and ML 
problems are nonconvex and difficult to solve. To avoid difficulty in solving the ML (or LS), authors 
in iflOl took a geometric interpretation into account and formulated the positioning problem as a CFP 
and applied the well-known POCS approach to solve the positioning problem. 

To formulate POCS, note that in the absence of measurement errors, i.e., di = dj(x, a^), it is clear 
that the target, at unknown position x, can be found in the intersection of a number of spheres with radii 
dj(x, a.i) and centers For nonnegative measurement errors, we relax spheres to balls and deduce that 




(15) 



i=i 



1 1 

the target definitely lies inside the intersection of a number of balls. Let us define the (closed bounded) 
ball Bi centered at &j as 

Bi = {x£R n : ||x-a*|| 2 < di}, i = l,...,N. (16) 

It is then reasonable to define an estimate of x as a point in the intersection B (a closed bounded set) of 
the balls Bi (a feasible point) as 

N 

x e B = P| Bi. (17) 

1=1 

Therefore, the positioning problem can be rendered to the following convex feasibility problem (CFP): 

minimize 

xeR" 

subject to ||x — aj || < di, i = l,...,N. (18) 
To solve (fT8T ). we note that CFP can be reformulated by minimizing the following convex function 

/(x) =max{dist(x,Bi),...,dist(x,Sjv)}, (19) 

with dist(x, Bi) denoting the minimum distance between x and any point in set Bi. 

Using negative subgradient updating method lfT2l . ||29l , we can obtain a solution to ( fT9l) by 

x fc+i = x fc _ afcg fc ; fe = 0, 1, . . . , (20) 



where x is the kth iterate, at is the /cth step size, and g is a subgradienu. A subgradient g of / at 
x fc can be computed as 

0, if /(x fc ) = 0, 

(21) 

l| x X ^°ff)|| a , if /(x fe ) + 0, dist(x fe ,^-) > dist(x fc ,^), Vi + j, 
where Vs j ('x k ) is the orthogonal projection of x fc onto the set Bj. By choosing the step size as = 
/(x fc )/||g fc ||2 in d20l , according to Polyak approach |[T2l . we derive the following approach, called 
alternating projections ||30l or POCS, for updating 

x fc+1 =P Bj (x fe ), A; = 0,1,..., (22) 

2 Let P be a nonempty set in K". A vector g G 1" is a subgradient of a function / :D->latx€fif /(y) > 
/(x) + g T (y - x) for all y 6 V m 
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Fig. 3. A 2-dimensional network consisting of three reference nodes and one target node. For nonnegative measurement errors, 
the target node at position x is found in the intersection of three discs. The POCS estimate converges to a point x inside the 
intersection area (in this case on the boundary). 

where index j is the one used in (|2TT) . 

As mentioned before, POCS gives an estimate that is feasible (if the intersection B is nonempty). In 
each step, POCS projects the current point x fc onto the farthest convex set. For example, Fig.[3] shows 
a 2-dimensional network in which the measured distances in reference nodes are positively biased. The 
POCS' estimate in this figure converges to a point in the intersection of three discs after two iterations. 
For more details on variations of the POCS algorithm and the application of POCS for the positioning 
problem, we refer the reader to lfT2l and O, ifTTll . |[3ll . respectively. 

V. Geometric upper bounds 

In this study, taking the assumption of positively biased measurement errors into account and consid- 
ering discussions in Section ITl-CL we derive two different upper bounds. The first bound is derived based 
on the availability of an estimate. If such an estimate is available (feasible or infeasible), we can bound 
it by finding the maximum distance between the estimate and any point in the feasible set. The second 
bound is derived without the need for an estimate, as the maximum length of the intersection set. 
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Let us define the norm of position estimate, which we call the position error, as 



l x - x l|2, 



(23) 



where x is an estimate of the target node position given by a positioning algorithm. In a practical scenario 
it is not possible to compute the exact position error in d23l since the position of a target node is unknown. 
Therefore, we may compute a lower or an upper bound on the position error for evaluation of an estimate. 
According to discussions in Section ITl-C I it seems that the plausible definition for the maximum position 
error, when a single estimate is available, can be considered as 



e < Umax,l = max ||x - x|| 2 , 



(24) 



where B defines a set (closed bounded) in which the target node x belongs. In fact, definition (1241 is a 
special case of the upper bound defined in ([8]) in Section Hl-CI when Af(m) = B. In other words, (124b 
defines the largest distance from a point to a set. 

Alternatively, if an algorithm always produces one point in the feasible set B as an estimate, we are 
still able to define an upper bound on the position error, even without having access to an estimate, by 
setting X(m) = X(m) = B in <[T0]>, 



e < t>max,3 = max ||x - y|| 2 . 

x,yGB 



(25) 



A. A bound for the case an estimate exists 

As mentioned in previous section, we can upper bound the position error due to an estimate x (either 
feasible or infeasible), by solving the optimization problem (1241) . The solution is found on the boundary 
of set B. For example, let us consider Fig. [4] where an estimate x of the target node position inside the 
intersection of three discs is available. The position error and the maximum position error are shown in this 
figure. Instead of directly solving the problem in (1241 . we consider a QCQP problem QP{Aj, bj, Ci}f =0 , 
where 



— I n , bj — < 



-x, if i = 0, 



-a,-, otherwise, 



Ci = < 



Ixll 2 



if i = 0, 
a i\\ 2 ~ df-, otherwise. 



(26) 
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Fig. 4. The position error and the maximum position error for an estimate x of the target for the network considered in Fig.[3] 



Obviously, t> qp {Aj, hi, Ci}f =Q = x . The optimization problem in (l24l ) is nonconvex which makes the 
problem complicated. To solve the problem, we employ a relaxation technique. Following the procedures 
explained in Section Hl-B I we can get a relaxed SDP problem as SDP{Aj, bj, q}^ and the maximum 
position error can be upper bounded as 

e = ||x - x|| 2 < u maXi i < yj v sd p{Ai, hi, Ci}f =Q . (27) 

In order to investigate the tightness of the upper-bound derived in (127T ). we can derive a lower-bound on 
-Uq P {Aj, hi, Ci}jL . Let us write the QCQP problem QP{A,;, hi, Ci}f =0 parameterized in (l26l ) as 

maximize tr ( B Tx T t] T Tx t t] 

xgR", rGR V l j l 

subject to tr (b^ [x t t] T [x T r] ) <t h i = l,...,N+ 1, 

where 



(28) 



B 



iV+l 



ti = $ + e 2 , i<N, t N+1 = l, 



Onxn On 




In -X 




In a i 




, B = 




i Bj — 




T 1 

"n ■*■ 




-x T ||x|| 2 




— || a i ||| "1" g2 



(29) 
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where e 7^ is any nonzero real value. It is seen that Bi y for 1 < i < N. Then, B * >~ °> 

meaning the interior of the feasible set is nonempty. 

Proposition 5.1: A lower bound on the optimal value of QP{Aj, hi, Ci}fL parameterized in (1261 ) based 
on the optimal value f s d p {Aj, b«, Ci}fL , can be obtained as 



/a^sd P {Ai,bj,Ci}^ < v qp {Ai,hi,Ci} i=0 , (30) 
where 

a= 21n(2(^ + lV) ' M = min{iV + l,n + l}. (31) 



Proof: Recalling the results of 11321 , which determines a lower bound on the optimal value of a 
QCQP based on its relaxed SDP, we get a lower bound on the optimal value of d28l ), which is exactly 

Vqp{A-i, bj, Cj}^_g, as 

aw s dp{Ai,bj,Q}^ < Vq P {Aj,bi,Ci}^ , (32) 

where 

a= 2H2(N + lW M = ^W^ + l,m f crank(B i )}. 
It is clear that rank(Bj) = n + 1. Therefore, a lower bound on f qp {Aj, bj, Ci}^L Q can be derived as ( |30l ). 

■ 

For details of deriving lower bounds on a nonconvex QCQP, we refer the reader to |[T8l . |24ll . |[32l and 
references therein. 

B. Bound regarding the feasible set 

In this section, we investigate another upper bound defined in (T25T ) and repeated here for convenience 

«max,3 = max{||x - v|| 2 : x,v€#}. (33) 
If a feasible point x G B is available, it is expected that the first upper bound w max ,i yields a tighter 



bound compared to the bound defined in (1331 ) (the maximum length of the intersection). In fact for a 
fixed xe6, 

max llx — wlU > max llx — xlU. (34) 

x,wGB xGS 
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Fig. 5. Maximum Euclidian distance of the intersection as the maximum position error for an estimate inside the intersection 
area of Fig.[3] 

The optimization problem in (l33l is nonconvex. Geometrically, it can be imagined as the diameter of 
the minimum ball enclosing the intersection. Instead of solving the problem formulated in (l33T ). we find 
a minimum ball covering the intersection B. Let us consider the center x c and the radius R of such a 
ball and formulate the minimum ball enclosing the intersection B in decision variables x c and 7 = R 2 
as 

minimize 7 

x c £M", 7SM+ 

subject to ||x-x c || 2 < 7, x 6 B. (35) 

Let the optimal solution of d33T ) be v' max 3 . Then, f max ,3 = 2y / t> max 3- Fixing x c in (l35T ). using Lemma |2T] 
and following a similar approach as used in |[33l . we can get the following optimization problem to find 
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the minimum ball enclosing the intersection B: 



minimize 7 

7GK+, AGR™ 



subject to 



(Sill \. ~ X c — YliLl ^i^i 



to, 



(xc - £f=i A ia ,) T 7 - ||y||5 + E-Ii MINI! - < 2 ) 

x c G K n . (36) 

Taking similar steps as done in |[33l . which implies for the optimal solution Yli=i Aj = 1 and x c = 
Si=i Ajaj, we can obtain an optimization problem to find an upper bound on the squared radius of the 
minimum ball enclosing the set B in the Euclidian norm sense as 

N N 

minimize || V" - Ai(||ai||| - 

T 2=1 1 = 1 

N 

subject to Aj = 1. (37) 

i=l 

Finally, an upper bound on the maximum length of B is given by 

%ax,3 < 2i?, (38) 



where R = yj\\ £^ A,a,||l - E,=i A*C[|«*II3 " <?)■ 

It has been proved in ll33l that when the number of constraints TV (here the number of reference nodes) 
is equal or less than n (the size of dimension), (l37l) gives the optimal solution to (l35l ). Otherwise when 
N > n, the optimal solution in (l37l) is an upper-bound to the optimal solution in (I35T ). The upper bound 
obtained by solving (l37l) then gives the maximum Euclidian length of the intersection. 

Another approach to compute an upper bound on w max ,3 is to replace B with an enclosing set in (|25T ). 
We will in the following consider two such sets. The first enclosing set is the bounding box|j for B, and, 
given the bounding box, it is very easy to compute an upper bound on t> ma x,3> see Fig. [6] 

The second enclosing set is found be replacing B- L with their bounding boxes, i.e., the £2 balls in (fT6l ) 
are replaced by the corresponding balls, 

Bi' = {xeR n : Hx-a.Joo <d t }, 

3 By the bounding box of the set A, we mean the smallest cuboid 1341 that is enclosing .4. 
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Fig. 6. The maximum length of the bounding box of the intersection as an upper bound for the network considered in Fig.[3] 



and noting that 

N 



B c B' = p| Bi'. 



i=i 

Hence, an upper bound to u ma x,3 is found by considering the length of B', see Fig. [7] 
To compute the bounding box for B, we study the following optimization problem: 

maximize ||x — y||oo 

subject to x, y G B. (39) 
The optimization problem in d39l again is nonconvex. Using the definition of the norm, we can write 

maximize max(|xi — y±\, . . . , \x n — y n \) 

subject to x,y e B. (40) 
The max function in d40l ) can be computed as 

max{ai, . . . , a n } = ai <s=^> a, > aj, Vj. (41) 
Using a dummy variable (3, we have 

maxjcti, . . . , a n } > /? <^=^> a± > /? or ai > (3 ... or a n > /3. (42) 
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Thus, using a simple technique, we need to solve two optimization problems for every dimension I as 
follows: 

maximize B 

/3GR 

subject to ||x — aj|| < dj, i = 1,... ,N, 

x e > 8, (43a) 

minimize B 

/3eR 

subject to ||x — aj|| < di, i = 1, . . . , N 

x e < B. (43b) 



The optimization problems in (1431 are called the second order cone program which is a special case 
of the quadratic programming. It can be easily transformed to an SDP lfl7l . Suppose that the optimal 
solution to problems (I43ab and (I43bb along a dimension I are x\ and , respectively. Let the maximum 
length for the fth dimension be t> SO cp,£ = \ x *i 1 ~ x e 2 \- Then, the maximum length of the intersection can 
be upper bounded as 



, ^socp,^) 2 . ( 44 ) 
\ i=i 



Thus 

^max,3 — ^socp- (45) 

To compute the upper bound on v maXj 3 based on B' , we consider the following optimization problem: 

maximize ||x — y||oo 
x,y 

subject to x, y £ B , (46) 

For example Fig.|7] shows the concept of relaxing the constraint for a 2-dimensional network. Following 
the same procedure to obtain (l43l . we obtain two optimization problems, called linear programs (LPs), 
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Reference node 
Target node 




Fig. 7. Every constraint is replaced with a bounding box and then a bounding box enclosing the intersection of relaxed 
constraints is computed. The maximum length of the bounding box enclosing the intersection gives an upper bound for the 
network considered in Fig.[3] 



for every dimension. For instance, two LPs for the £th dimension can be written as 



maximize 



subject to ti — cii : e — di < 0, 

t e - au + ck < 0, i = l,...,N, 



(47a) 



minimize 



subject to ti — an — di < 0, 



<H4-te + di<0, i = l,...,N. (47b) 
The optimal solution to the optimization problem (|47T ). i.e., t% and t\ , are simply computed as 

t* £i = min{ai i £ + d\, . . . , a N j + d N }, t} 2 = max{ai^ - di, . . . , cln,i - div}. (48) 



Let vipj = \t* i± - t\ 



1, . . . , n, be the maximum length along the £th dimension. The maximum 
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TABLE I 

Summary of bounds. 



Definition: 


Eqn. 




e = x - x 2 


d23} 




Dmax,l = max X - X 2 


$24$ 




«max,3 — max x - y 2 

x,y£B 


{25} 




Upper Bounds: 




Eqn. 


Bound 1: 






e < «max,l < \/v s d p {Ai 


h r\ N 




Bound2: 






"max, 3 < 2R 




138} 


Bound3 (Type 1): 






Vmax,3 ^ Vsocp 




(|45) 


Bound3 (Type 2): 






■Umax, 3 < «lp 




us 



length of the intersection B is then upper bounded by 



V 



Em- 



(49) 



Therefore, an upper bound on position error based on a bounding box approach is given by 

Wmax,3 < (50) 

It is clear that f SOC p < %v 

TableJ] summarizes the various types of bounds derived in this study. 

VI. Simulation results 

In this section we evaluate the validity of different upper bounds. We consider a 1000 m 3 cubic space 
for simulation. N reference nodes are randomly distributed in the space. One target node is randomly 
placed inside the volume. To add measurement noise to actual distances between reference and target 
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nodes, we use an exponential distribution denned as 

/(*) = < 

0, ej < 0. 

The mean 1/7 is set to 1 m. The validity of exponential distribution, especially for NLOS conditions, has 
been justified in the literature, e.g., ifTTI . 1261 , ||35l . We study the POCS algorithm that always gives an 
estimate inside the intersection B in ( fT71 ) . To solve the optimization problems formulated in this study, 
we use the CVX toolbox 1361 . 

To evaluate the tightness of the bounds in Table H we consider the normalized difference between a 
bound v and the true error e, i.e., (v — e)/e. To illustrate how the tightness varies with, e.g., network 
deployment, measurement noise, estimator parameters, we study the cumulative distribution function 
(CDF) 

P v (x) =Pr j^-^ < x 

where the randomness comes from selecting, e.g., the deployment in a random fashion. In the following, 
we will generate e from POCS estimates. Since an estimate of the target position is available, we also 
consider the first upper bound for further comparisons. In all simulations, we generate 1000 random 
networks. 

Fig-El shows the CDF of the normalized position error of an upper bound versus POCS position error 
for different number of reference nodes. As expected, the first upper bound shows better performance 



compared to the other bounds. For instance, Fig. 8(a) shows that in 80% of the cases, the first upper bound 



computed by the network consisting of five reference nodes is less than 2.3 times the actual position 
error (considering the normalized error (v — e)/e). This figure also shows that the upper bound 3 (Type 
2) is the loosest bound. When the number of reference nodes increases, the upper bound 3 (Type 1) gets 
closer to the upper bound 2. Roughly speaking except for the upper bound 3 (Type 1), we can say that the 
behavior of other upper bounds (based on the normalized error (v — e)/e) does not change considerably 
with increasing the number of reference nodes. Fig.[8] also shows that the proposed bounds always are 
upper bounds (although not always tight). 
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In the next simulation, we compare the upper bounds with the maximum position error. To compare 
four upper bounds, we again employ the POCS method. For every realization of the network, we run 
POCS for 200 random initializations and take the maximum position error. For every realization, the 
upper bound 1 corresponds to the maximum distance to the intersection for the estimate that gives the 
maximum POCS position error. Three other bounds are independent of the POCS estimate and they 
approximate the maximum length of the intersection area for every realization. Fig. [9] plots the four upper 



bounds against the maximum POCS position error. In Fig. |9(a)[ we plot the upper bound 1 and a lower 
bound on the maximum position error when an estimate is available. As seen, the maximum position 
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POCS maximum position error [m] 




POCS maximum position error [m] 
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(b) 



4 6 8 10 

POCS maximum error [m] 



(c) (d) 

Fig. 9. Comparison between three upper bounds and the maximum position error of POCS for 15 reference nodes and 200 
random initializations for every realization, |(a)| Boundl is computed using the estimate that gives the maximum position error 
for POCS, [(b)! Bound2, |(c)l Bound3 (Typel), and[(d)lBound3 (Type2). 



error is bounded between the green and black curves, which defines an upper and a lower bounds on the 
maximum position error, respectively. These figures graphically show that the upper bound 1 is tighter 
than other bounds. They also show that the upper bound 3 (Type 2) is the loosest one. 

In Fig.flOl we plot the CDF of the normalized position error of upper bounds versus the maximum 
POCS position error for different number of reference nodes. Roughly speaking, in more than 90% of 
cases the upper bound 1 is equal or less than 1.5 times the maximum POCS position error for different 
number of reference nodes. Again, we see that the upper bound 1 is the tightest and the upper bound 3 
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Normalized error Normalized error 

(c) (d) 

Fig. 10. Comparison between the CDF of normalized error of different upper bounds versus the maximum position error of 
POCS for, |(a)| 5 reference nodes, [(b)] 10 reference nodes, [(c)] 15 reference nodes, and |(d)| 20 reference nodes. 



(Type 2) is the loosest one. It is seen that when the number of reference nodes increases to 15, the upper 
bound 2 in 80% of cases is tighter that the upper bound 3 (Type 1). 

VII. Conclusions 

In this paper we have formulated a number of upper bounds on the realization of the positioning error, 
i.e., the error which is produced by an estimator, or a class of estimators, given a certain realization of 
the measurement, m. The bound defined in ([8]) can be computed by finding the largest distance between 
a point in the set X(m), i.e., the set of all possible positions of the unknown node, conditioned on the 
observation m, and the estimate x(m, f ). (Recall that f contains the estimation algorithm parameters, e.g., 
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initialization, that determines how m is mapped to the position estimate.) Similarly, the bound in (fTOl) can 
be computed as the largest distance between a point in X(m) and a point in X(m) = {x(m, f) : f € J 7 }, 
i.e., the set of all possible estimates in the class of estimators defined by T . Hence, the bounds are 
nontrivial (i.e., finite) only if the measurement implies that the above-mentioned sets are of finite lengths. 
Moreover, it is, in general, not clear if the bounds can be computed with reasonable complexity. 

However, we have showed that we can indeed compute nontrivial bounds in an efficient manner for 
the special, but interesting, case when m consists of positively biased distances estimates between a 
number of reference (anchors nodes) at a-priori known positions and a target node (at an unknown 
position). We note that non-negative distance errors are likely to occur in non-line-of-sight environments. 
For this special case, the target node is constrained to be in the intersection B of a number of balls, Bi, 
i = 1,2, . . . , N, which are centered around the reference nodes and whose radii are given by the observed 
distance estimates. That is, in this special case, X(m) = B. An efficient algorithm, (|27T ). can then be 
found by relaxing the original bound d24l ) into a convex optimization problem using SDP techniques. 

Moreover, if we use a POCS algorithm to estimate the target node position, we know that X (m) = B, 
i.e., the estimate will be in B. Hence, the bound ((8]) simplifies to d25l ). To arrive at bounds that can be 
efficiently computed, we formulate three upper bounds of (|25T ) in (l38l ). (l45l ). and (l50l ). The bound (I38T ) is 
based on SDP relaxation, the bound (|45T ) by replacing B with its bounding box in (125T ). and the bound (l50l ) 
by replacing Bi with their bounding boxes in (fTTT ). Simulation results based on the POCS estimate for 
different situations show that the proposed upper bounds provide reasonably tight bounds. As expected 
from the theoretical part and confirmed by the simulation results, for the POCS estimate the first bound 
in (|2VT > is the tightest bound among different upper bounds formulated in this paper. The numerical results 
also show that the behavior of different bounds, except the one in d45l ), based on the normalized error 
does not considerably change with node density. It is also concluded from both theoretical aspects and 
simulation results that the bounds (|38l ) and (|45l ) are tighter than the one in (|50l ). 

Finally, it is clear that it is very valuable if we, in a practical situation, can append an estimated 
position with an upper bound of the position error. This is much stronger than saying something about 
the statistics of the position error (e.g., the mean squared error). The methods developed in this paper 
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provides tools for bounding the position error, albeit in somewhat limited situations, i.e., when X(m) 
has finite length. There are practical situations where this is a valid assumption, but also cases when it 
is not. 
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